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We describe a model of jet quenching in nuclear collisions at RHIC energies. In 
the model, jet quenching is to be caused by the interruption of jet formation by 
nucleons arriving at the position of jet formation in a time shorter than the jet 
formation time. Our mechanism predicts suppression of high-pT spectra also in 
d+Au reactions. 



PACS: 25.75.-q, 25.45.-z, 25.75.Nq 

1 Introduction 

Hadronic spectra at high transverse momenta from uhrarelativistic nuclear colhsions re- 
sult from jets produced in the early hard partonic interactions. The yields are normally 
expected to scale with the number of binary nucleon-nucleon collisions. However, ex- 
perimental results obtained at the Relativistic Heavy Ion Collider (RHIC) contradict 
these expectations: in central Au-f Au collisions the yields above 2 GeV/c are strongly 
suppressed with respect to binary collisions scahng Pl^lEl^. Recent data from d-|-Au 
collisions show an enhancement of the yield with respect to the binary-collision scaling 
for pt ^2, GeV/c [3 El El but indicate that the yield may fall below this expectation 
above 8 GeV/c 0. 

The lack of strong suppression of high-p^ yields in deuteron-induced reactions sug- 
gests that the effect in Au-f Au collisions is due to partonic energy loss in a medium 
with very high energy density. This mechanism has been studied by many authors 
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In addition, there are some models which aim to describe high-p^ spectra from nuclear 
collisions by invoking different physical effects. The perturbative-QCD-improved parton 
model of refs. jElE] includes shadowing and the broadening of spectra due to intrinsic 
parton transverse momentum. The authors conclude that the mechanism is insufficient 
to explain high-p^ suppression at RHIC and even at the SPS. 

The mechanism of parton saturation ^1 ^1 was argued to lead to suppression of 
intermediate pt production in both Au+Au and d+Au collision systems POII^ . 

In ref. |22| . we proposed a model for the suppression of high-pj^ spectra in ultrarel- 
ativistic nuclear collisions based on arguments involving the uncertainty principle and 
formation time [2l^\ I24j . The mean-free path of an incident nucleon in the nucleon- 
nucleon centre-of-mass system (CMS) at RHIC is rather short: A « 0.025 fm. When two 
jets are produced in a particular nucleon-nucleon collision, the next nucleon arrives at 
the position where the jets were created in a time interval of the order of A/c. According 
to the uncertainty relation, a process with longitudinal momentum transfer Ap^ and 
energy transfer AE needs space and time intervals given by 

h , h 

> -r , At > — . 1 

ApL AE ^ ^ 

to reach completion. Hence, if the mean-free path or the mean-free time are shorter 
than these intervals, the created jets may be seriously influenced by nucleons arriving at 
the position where the process develops. For the quoted mean-free path at RHIC, this 
puts a limit on processes which are not influenced by our mechanism at ApL > 8 GeV/c. 
The limit is, of course, not strict, since the time interval between subsequent collisions 
fluctuates. 

In the next Section we present a slightly improved version of the model from our 
earlier paper [22 . The results are shown in Section O and we conclude in Section^ The 
appendix contains a discussion of Lorentz invariancc in our model. In what follows we 
shall work in natural units c = fi. = 1. 



2 The model 

We shall formulate our Glauber model in the nucleon-nucleon centre-of-mass system. For 
the sake of simplicity, we assume that the density distribution is uniform. Technically, 
the model is rather similar to those for nuclear absorption of J/ip in heavy-ion collisions 

m 

A Glauber model is usually formulated in terms of "tube-on-tube" collisions, with 
tubes filled with nucleons^. For a collision of nuclei A and B at impact parameter b we 
have for the lengths of the colliding tubes 

2La(s) = 27-iy^i?^ - s2 , 2LB{b,s,e) ^ ^ Rl - b'^ - + 2b scose . (2) 

^The integral of the density along the "tube", i.e. in the longitudinal direction, is often called the 
thickness function = p{%/s^~+~z^) dz. In our model, the density is given by paS(Ra ~ '") a^nd 

the thickness function corresponds to the tube length given in eq. |2j multiplied by the constant nuclear 
density. 
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Fig. 1. Left: geometry of non-central nuclear collisions. Right: layout of tube-on-tube interac- 
tion (plotted without Lorentz contraction). 

Here, 7 is the Lorentz contraction factor for the boost of the nuclei from their own rest- 
frame to the CMS. All other coordinates and sizes on the r.h.s., however, are taken in 
the rest frames of the nuclei and their meaning is explained in Figure ^ The positions 
of nucleons within the colliding tubes will be denoted by za and zb- The values of za 
and Zb satisfy 

-La < za < La, -Lb < zb < Lb- 
We take za and as increasing in the direction of motion of A and B respectively in 



We will be interested in comparing the yield of jets at high transverse momentum 
produced in nuclear collisions to that produced in nucleon-nucleon collisions. By staying 
at the level of jets and not including fragmentation, our calculation is not influenced by 
the form of the fragmentation function which is not entirely known for nuclear collisions. 
On the other hand, we must keep in mind that our results cannot be compared directly 
to the measured pT-spectra. 

The yield of jets in collisions of nuclei A-|-B at impact parameter b is defined as 



the CMS. 



Yab{ptM 



(3) 



dcTAB 

db^ 



while the yield in nucleon-nucleon collisions is introduced as 




(4) 



cr, 



PP 
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As in our previous paper j22| we intend to determine quantity 



YppiPT) 

p pLa pLb 

Gnn / sdsdO I dzARA / dzB PB F {b, S , 9 , ZA, Zb) ■ (5) 

^overlap J —La — Lb 

Here, pA and pB are the Lorentz-contracted nuclear densities of the nuclei A and B, 

PA = PB = P = 1 Po 

with po = 0.138 fm^"^. For the total inelastic nucleon-nucleon cross-section we will take 
o'nn ~ 40 mb = 4 fm^ . 

Without the suppression factor F{b, s,6, za, zb), formula lO would give the average 
number of binary collisions in an interaction of nuclei A+B at impact parameter b. The 
suppression factor will account for the assumed effect of jet destruction. 

If jets are created in a hard interaction of two incident partons of the colliding nuclei, 
we shall assume a cross-section for their destruction by a subsequent incoming nucleon 
(in CMS) an expression 

"°^^^'^'^^^^"°( i + (Mt--t))0 ' 

where t is the time coordinate of the hard interaction and t' is the time^ when the 
destroying nucleon arrives at the position of the jet"*. The jet transverse momentum px 
is roughly equal to the energy AE involved in the process. Therefore, the form is in 
line with our considerations about the jet destruction based on the uncertainty relation 
described in the introduction. The absorption cross-section is of order of few milibarns 
and were tuned by the parameter erg . Although formula JHl for the cross-section is written 
in the CMS reference frame, we demonstrate in the Appendix how it can be written in 
an explicitly invariant way. The improvement over our previous paper |22| lies in the 
time-dependent prescription for the absorption cross-section. In ,5^ we ignored the fact 
that nuclcons with different distances from the origin of the jet have different chances to 
destroy it, and we assumed 

a^ipr) = ao ( ^ ^ (p^'/^^^). ) ^ (^^^ -odel of M) , (7) 



^In our previous paper 1221 we used the longitudinal coordinates z and z' instead of the times in our 
argumentation. The advantage of the new formalism is that cTq can be written in a Lorentz invariant 
way, as described in the Appendix. At the energies studied, however, the numerical difference between 
{t' — t) and (z' — z) is negligible. 

*For brevity, we talk about a "destroying nucleon" although we rather mean those of its quarks with 
large enough momentum. The problem with soft sea quarks is that they can be barely localised as 
precisely as we need for our argument, but we can assume that their energy is too small to destroy the 
jet. 
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with ptq = 8 GeV/c for collisions at RHIC. 

The suppression factor F{b, s, 6, Zj^, zb) now reads 



s, 6, ZA, Zb) = exp ^- ^ Oa (pt, t, ^) pa v dtj 

(- [ aa(pT,t,—] pBvdt] 

\ J-Ls/v ^ V ^ J 



X exp 



(8) 



where za and zb determine the position of the hard interaction within nuclei A and B, 
respectively, and v is the longitudinal velocity of the colliding nuclei. The b, s, and 0- 
dependences of F are implicitly included in La and Lb. When we insert the prescription 
(01 into equation the integrals can be performed analytically and lead to 



F{b,s,0,ZA,ZB) = 

<^ovpA ( J A 

exp 



r, I 1 , +arctan/^ 



2pT 



1 + 



arctan /b 



with 



Ja ^—{za + La) , 

V 



fs = —{zb + Lb) ■ 

V 



(9) 



(10) 



We can compare our results with our previous model of ref. (22j . The suppression 
factor there was obtained by inserting formula for the absorption cross-section. This 
leads to 

F°'''{b,s,0,ZA,ZB) - cxp(^-J'^ a:'''{pT)pAdz^exp(^- J af''{pT)pBdz^ 

= exp {-a°"'{PT) PA (zA + La) - a^'^ipr) Pb (zb + Lb)) . 

(11) 



3 Results 



In our calculation we simulated a set of nuclear collisions within the Glauber model 
framework and for every nucleon-nucleon interaction we determined the suppression 
factor according to eq. By adding the suppression factors from all nucleon-nucleon 
interactions we obtained the total yield of jets relative to the yield from nucleon-nucleon 
interactions as defined in eq. The set of nuclear collisions under study was always 
chosen appropriately to sample given centrality requirements. 
The absorption cross-section parameter ctq was tuned to 

do = 8mb = 0.8 fm^ . 

In the figures we always plot the relative yield of jets divided by the number of 
binary collisions. In case of no suppression this number approaches unity. As we have 
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Fig. 2. The ratio -Rau+Au(pt, &)/?^coll(&) for central collisions at ^/s = 200 ^IGeV plotted as a 
function of ncoii- Different curves correspond to variations of pT, as indicated. 



already mentioned, we cannot compare our results directly to experimental data because 
fragmentation is not included in our calculation. Nevertheless, we make some remarks 
on the relation of our results to data in Sectional 

Figure|2shows the dependence of the relative yield on the number of binary collisions. 
We observe that scaling with the number of binary collisions is slowly recovered for large 
Pt- As the transverse momentum increases, the curves converge to the asymptotic value 

of i?AB(PT, "coll )/"coU = 1- 

The pT-dependence of the yields of jets is plotted in Figure|21for central and peripheral 
Au+Au collisions and for d+Au collisions. Our mechanism leads to suppression of jet 
production even in d+Au collisions. We comment on the relation of this result to data 
in the next Section. 

In Figure 01 we compare our model with that of ref. |22| . This comparison shows 
that the models give similar results at low transverse momenta while at high "pj- the 
suppression is stronger in the new model. This is due to destroying nucleons which may 
be very close to the jet production site in the new model (i.e., (i! — t) is small) and thus 
associated with a very large absorption cross-section. In the old model , on the other 
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Fig. 3. The ratio i?Au+Au(pT, 6)/ncoii(t') plotted as a function of the jet transverse momentum 
Pt- Calculated for collisions at ^/s — 200 AGeV. Different curves refer to: central Au+Au 
collisions (0-5% of the total cross-section), peripheral Au+Au collisions (60-80%), and minimum 
bias d+Au collisions. 



hand, all nucleons were effectively put a distance A from the jet production such that 
there were effectively no nucleons that would be close enough to destroy very-high-pr 
jets^. In Figure 0] we also show the result of a modified model in which we forbid the 
nucleons within a nucleus (at rest) to be closer to each other than 1 fm, such that there 
are no very close destroying nucleons. This makes the suppression weaker by 5-15%, 
depending on pt- 

When we compare the old and the actual model at low pt, we see the effect of 
suppression by subsequent nucleons folowing the first potentially destroying nucleon. 
While in the old model their absorption cross-sections were equal to that of the first 
destroying nucleon, in the new formulation these cross-sections are smaller due to the 
larger distance from the origin of the jet production. Therefore, the new model gives 



^Mathematically, in the present model the supression factor F in eq. j5J goes like exp(— const /p|,) 
for high Pt, while in the old model the factor F°^'^ of eq. 1111 approaches exp(— const /p^). 
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Fig. 4. The comparison of our model to that of our earlier paper |22|. The relative yield 
RAu+Au{pT,b)/ncou(b) plotted as a function of pt, calculated for central Au+Au collisions (0- 
5% centrality cut) a.t = 200 ylGeV. 



weaker suppression at low px than the old one. 

4 Conclusions 

We have described here an improved version of the jet attenuation model of |21j in which 
the cross-section for the jet attenuation is time-dependent as implied by We have 
shown that results for the jet attenuation obtained by the new version do not differ very 
much from the previous version. 

An important feature of our model is the lack of any thresholds under which the 
mechanism of jet suppression would cease. We obtain a suppression with respect to the 
scaling with the number of binary collisions even in the case of d+Au collisions. This 
differs from the recently published high-pT data from d+Au collisions at RHIC [S1|S1[7||SI- 
An enhancement, i.e., a yield larger than expected from the ncoU-scaling was reported 
from the deuteron-induced reactions. This is likely to be, at least in part, a Cronin effect 
which results from pT-broadening of the incoming partons. Since there is no such effect 
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included in the model, we cannot describe this feature. 

However, the data on high-pT spectra of charged hadrons from d+Au collisions pub- 
lished by STAR collaboration indicate [H] that the yield normalised to the pr-spectrum 
and the number of binary collisions possibly falls below unity for pt ^ 9 GeV/c. In this 
kinematic region, the Cronin effect may not be effective any more and the suppression 
could be described by our model. In order to indicate this, in Figure |21 we did not plot 
the curves in the kinematic region where we expect the dominance of the Cronin effect 
and highlighted only the kinematic region where our model could stay relevant as is. We 
plan to combine the Cronin effect with the present model in the future and confront our 
results with existing data. 



One of the authors (JP) is indebted to the CERN Theory Division for the hospitality 
extended to him. We like to thank Orlando Villalobos-Baillie for many useful discussions 
and careful reading of the paper. The work of JP and NP was supported in part by the 
Slovak Ministry of Education under the Grant No. VEGA V2F13. 

A Lorentz invariant formulation of the absorption cross-section 

In the earlier version of our model |23, the relative yield R, the suppression factor 
F, and the absorption cross-section Ua, given by eqs. (6), (7), and (3) of that paper, 
respectively, are manifestly invariant under longitudinal boosts, since pr and (Tq are boost 
invariant and the products pA dzA and ps dzB are as well. In the present formulation, 
the corresponding equations (0), JHl, and ^ are valid in the CMS frame, but they can 
be also written in a Lorentz invariant way. 

Prescription © for the absorption cross-section was formulated in the centre-of-mass 
system. In this system, it says that (Ja will be large if 



We treat only jets at central rapidity which do not move longitudinally in CMS. The 
positions where the jet is created and where it is hit by another nucleon are thus the 
same. Thus, the l.h.s. of the previous inequality can be written in a Lorentz invariant 
way: 
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(12) 




(13) 



where x'^ and x"^ are the four- vectors corresponding to events of jet creation and possible 
destruction, respectively. A Lorentz invariant formula for the absorption cross-section 
for mid-rapidity jets thus reads 
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We can put this in a more suitable form. We neglect the movement in transverse 

direction and write 

(a:"'-x'')(:c;,-x^)=Af2-Az2 (15) 

We denote: 
Uj the jet rapidity in a given frame, 

yA the nucleon rapidity, same as nucleus rapidity, 

z* the distance of the destroying nucleon from the place where the 
jet was produced, measured in the rest-frame of the nucleus. 

Then, the time between jet creation and its possible destruction can be written as 

Ai = — , -, ^ , — r , (16) 

cosh?/^(tanhyA — tanhj/j) 

and the distance between these two events is given by 

z* tanhvi , „, 

Az = At Vi^t = —r 7-^ \ , , • 17 

coshj/A(tanhyA — tanh?/j) 

Using these equations and after some algebra, we obtain 

At^ - = _J . (18) 

smh {yA - yj) 

This form depends only on difference of rapidities and is thus Lorentz invariant. 
If we are only interested in midrapidity jets, we can replace 

Va + Ub 

where yA and ys are rapidities of the colliding nuclei and obtain the absorption cross 
section 



1 + sinh^((yA-?;B)/2) , 
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